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We report the investigation of the generation and detection of GHz coherent acoustic phonons in 
plasmonic gold nanoparticles superlattices (NPS). The experiments have been performed from an 
optical femtosecond pump-probe scheme across the optical plasmon resonance of the superlattice. 

Our experiments allow to estimate the collective elastic response (sound velocity) of the NPS as 
well as an estimate of the nano-contact elastic stiffness. It appears that the light-induced coherent 
acoustic phonon pulse has a typical in-depth spatial extension of about 45 nm which is roughly 4 
times the optical skin depth in gold. The modeling of the transient optical reflectivity indicates 
that the mechanism of phonon generation is achieved through ultrafast heating of the NPS assisted 
by light excitation of the volume plasmon. These results demonstrate how it is possible to map the 
photon-electron-phonon interaction in subwavelength nanostructures. 


PACS numbers: 78.67.Pt, 73.20.Mf, 78.47.D-, 62.25.-g 

Plasmon assisted subwavelength light transmission 
is a key physical mechanism for modern nano-optics 
and nano-plasmonics [I.. The propagation of plasmon- 
polariton waves in nano particles nanostructures (chains, 
arrays) is at the core of intense fundamental investi¬ 
gations and has led to numerous reports [2®. For 
plasmonic applications, it is obviously crucial to un¬ 
derstand and control the damping of these plasmon- 
polariton waves, i.e. the collective propagation distance. 
This propagation distance as well as a more general de¬ 
scription of the dispersion curves for both longitudinal 
and transverse light electric field have already been ob¬ 
tained for supported ID chains |2H5|. The propagation 
distance is currently limited by the intrinsic electron- 
phonon, electron-defect interaction within each particle 
as well as radiation loss. In some particular 2D optical 
spectroscopy imaging, it has been shown recently that 
the propagation/distribution of the plasmon-polariton 
could even be mapped for ID nanoparticles chains 016]. 
The visualization of these plasmons at the interface 
air/nanostructure with near-field imaging [5] or with an 
electron beam imaging |7] were also reported. This state- 
of-art 2D imaging is however limited to ID or 2D systems, 
i.e. to surface plasmonics and no local measurement of 
volume plasmon-polariton propagation has been reported 
for 3D nanoparticles arrays so far. Most of the plas¬ 
monic responses in nanoparticles superlattices are ob¬ 
tained indeed from far field optical absorbance measure¬ 
ments nng. Because of the difficulty to probe the innner 
part of a plasmonic superlattices, no quantitative depth 
profiling description of the plasmon-polariton propaga¬ 
tion have been reported to date. However, it is known 
that ultrafast optical techniques can provide a picture 
in time and space of the light-matter interaction QIB- 
[12]. In particular, in case of femtosecond light pulses, 


the light-matter coupling can lead to the emission of co¬ 
herent acoustic phonons that occurs all over the spatial 
extension where the light energy is converted into me¬ 
chanical energy. Consequently, the characteristic spatial 
profile of the emitted coherent acoustic phonons contains 
information on the light-matter energy conversion pro¬ 
cess. In bulk solids, the analysis of coherent acoustic 
phonons packet has permitted to understand the evolu¬ 
tion in time and in space of the different electron-acoustic 
phonons interactions like the thermoelastic process (ul¬ 
trafast heating of the lattice) [T3MT5] . the deformation 
potential process (electronic pressure) [T6HT8] and even, 
reported recently, the spatial electron-hole separation in 
semiconductors (i.e. the photoinduced Dember electric 
field) jT9]. 

In this Letter, by applying optical pump-probe scheme 
depicted in Fig. [l] we report the observation of GHz 
coherent acoustic phonons generation in NPS. The laser 
excited coherent acoustic pulse has a spatial extension 
in the range of 45 nm which is much larger than either 
the gold nanoparticle diameter (3.7 nm) or the light pen¬ 
etration depth in gold. Because laser-excited hot elec¬ 
trons are confined in the gold nanoparticles, these obser¬ 
vations indicate that the acoustic phonons are very likely 
induced by ultrafast heating of the superlattice assisted 
by a collective plasmon-polariton propagation. While in¬ 
dividual vibrations of nanoparticles assemblies have been 
reported earlier mm or narrowband Brillouin mode 
in semitransparent cobalt superlattice [22], our results 
demonstrate that it is possible to generate propagating 
broadband coherent acoustic phonons in plasmonic NPS. 

The samples have been grown by soft chemical routes 
(for ample details see the review [23]). These nanopar¬ 
ticles, encapsulated with 0.7 nm long molecular chains 
(MSA Mercaptosucinic), are packed on a silicon sub- 
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Figure 1: (color online) (a) Transmission electron diffraction (TED) pattern showing the Debye Scherrer rings (top figure) 
coming from the crystallized cubic gold nanoparticles, tn the middle figure, the transmission electron microscopy image (TEM) 
shows the mesoscopic hexagonal arrangement (hep) with the two neighbor plane distance of 4.4 nm. The hep arrangement 
is also well evidenced with small angle electron diffraction revealing the 6-fold axis (bottom figure), (b) Sketch of the gold 
nanoparticles superlattices connected via MSA molecules, (c) Optical microscopy view of a superlattice deposited onto a 
silicon substrate, (d) AFM profile of the hexagonal shape superlattice (height H=180 nm). (e) Sketch of the pump and probe 
experiments performed on these nanoparticles superlattice where coherent acoustic phonons are generated and detected in the 
front-front configuration, (f) Top figure : plasmon resonance of the gold super lattice from [9] fitted with Eq. 0 . Middle and 
bottom figures : refractive index (real and imaginary parts) of the gold superlattice compared to those of bulk gold [27h The 
red symbols (square and circles) are the refractive index values estimated with the photoacoustic response (see Eq. pi and 
details in the text). 


strate according to an hexagonal structure (hcp-packing) 
with the c axis perpendicular to the substrate (the 6-fold 
symmetries are clearly visible in Fig. [lja)-(c)-(d) from 
mesoscopic to microscopic scales). The measured NPS 
spacing of d = 5.1 nm, gives the particle diameter of 3.7 
nm. Two different superlattices have been studied with 
typical lateral size of tens of micrometers and a thick¬ 
ness of 52 nm (triangular shape, Fig. [ljc)) and 180 nm 
(hexagonal shape, Fig. 0d)). These ordered assemblies 
exhibit plasmonic properties (Fig. 0f)) with a collective 
response modeled with the following simple phenomeno¬ 
logical dielectric function USES]: 


esi(w) = 1 + 


WZ 


; 2 — ^p/3 — cj 2 — iYuj 


(1) 


The real ( usl ) and imaginary (fcsx) parts of the re¬ 
fractive index, as well as the plasmon dielectric losses 
(oc 2nkuj [26] ) have been deduced from this dielectric 
function ( tisl = ns l + iksL = y/tsi) and is shown in 
Fig. [ljf) where us= 2.3 eV, cjp=2.2 eV and /rT=0.59 
eV, are the surface plasmon resonance energy of the in¬ 
dividual nanoparticle, the effective volume plasmon en¬ 
ergy and its damping energy respectively. One can notice 
that both the real and imaginary parts are clearly differ¬ 
ent from those of bulk gold crystal [27 (see middle and 
bottom panels in Fig. 0f))- By comparing this model 
and more recent numerical calculations [28], we can no¬ 
tice that Eq. 0 gives a reasonable approximation for 


very small nanoparticles (for instance 3.7 nm) while it 
does not for larger nanoparticles diameter (> 20 nm). 

The pump probe technique used here is based on a 
80 MHz repetition rate Ti:sapphire femtosecond laser. 
The beam is split with a polarizing beam splitter into a 
pump and a probe beams. The probe beam is introduced 
in a synchronously pumped OPO that allows to tune the 
wavelength and permit to perform two-color pump-probe 
experiments (pump and probe are linearly and circularly 
polarized respectively). The transient optical reflectivity 
signals are obtained thanks to a mechanical delay stage 
(delay line) which enables a controlled arrival time of the 
probe pulse regarding to the arrival of the pump pulse. 
The experiments were conducted with incident pump and 
probe beams perpendicular to the surface as shown in 
Fig. 0..)- 

The transient optical reflectivity signals obtained for 
superlattice are shown in Figs. [2] and [3] Time resolved 
reflectivity signals have been recorded for variable probe 
wavelengths (2.17-2.25 eV) but no significant variation 
have been observed in this wavelength range, probably 
because of the broad collective plasmon resonance (Fig. 
[TJf )). For the thicker NPS of 180 nm, the transient opti¬ 
cal reflectivity signal exhibits oscillatory components up 
to nearly 1 ns (pump 1.59 eV, probe 2.19 eV). In order 
to better reveal the oscillatory parts of the signal we pro¬ 
cessed the numerical time derivative shown in Fig. [2ja) . 
At first glance, the observation of the signal in the time 
domain reveals some periodic bursts repeated about ev- 
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Figure 2: (color online) (a) Time-derivative of transient op¬ 
tical reflectivity revealing the coherent acoustic phonon signal 
(blue curve) with numerical adjustment (red curve) contain¬ 
ing information on the acousto-plasmonic processes of the 
generation and detection of the acoustic phonons (see Eq. 
[ 2 ]) (inset shows the fast electronic response of the superlat¬ 
tice). As a comparison a simulation of the coherent acoustic 
phonons signal is given when the pump light energy is con¬ 
verted only within the gold crystal skin depth (i.e. 11 nm 
at 1.59 eV without hot electrons diffusion) and probed with 
standard optical properties of gold crystal (black curve), (b) 
Corresponding coherent acoustic phonons spectrum obtained 
by a Fast Fourier Transform (FFT) of the experimental time 
derivative transient signal shown in (a). Four harmonics of 
the mechanical resonance of the superlattice appear (inset). 



Figure 3: (color online) (a) Time-derivative of transient op¬ 
tical reflectivity revealing the coherent acoustic phonon signal 
(blue curve) with numerical adjustment with a damped sinus 
(red curve) (inset shows the fast electronic response of the 
superlattice), (b) Corresponding coherent acoustic phonons 
spectrum obtained by a Fast Fourier Transform (FFT) of the 
experimental time derivative transient signal shown in (a). 


ery 200 ps (Fig. [2|a)). We attribute this to an acoustic 
echo laser generated at the NPS free surface and trav¬ 
elling back and forth across the NPS and periodically 
detected at the front side upon reflection at the Si sub¬ 
strate interface. This interpretation is confirmed from 
the analysis of the FFT signals (Fig. lib)) where the 
detected frequencies follow the expected eigenmodes fre¬ 
quency sequence f n = tiVla/^H, with n=l, 2, 3... and 
Vla is the longitudinal sound velocity and H the NPS 
thickness (see inset of Fig. [2jb)). From the mean least 
square analysis of the slope, we extract a sound veloc¬ 
ity of Vla= 1770 m.s -1 which is about half of that mea¬ 
sured in bare gold crystals. The decrease of the sound 
velocity is an indication of the soft nanocontacts between 
nanoparticles. The time derivative signal shown in Fig. 
|3ja) corresponding to the thinnest NPS of 52 nm thick¬ 
ness resembles a damped sinus function of about 18-19 
GHz fitted frequency (Fig. [2](b)) . As already observed in 
different nanometric thin films [29H3T] . since the NPS is 
much thinner (52 nm) than the previous one (180 nm), 
the photoexcitation leads to the mechanical resonance of 
the entire superlattice layer at the first predominant har¬ 
monic frequency /1 which leads to the estimate of V la 
« 1870-1970 m.s -1 , in close agreement with the previous 
estimate obtained for the thicker superlattice. The esti¬ 
mate of the sound velocity provides a direct evaluation of 
the molecular contact elastic stiffness. Following the hep 
lattice dynamic equation along the z axis, we can connect 
the sound velocity to the nanoparticle contact effective 
stiffness K with V LA={^/V^)d^K/m [32] where d= 5.1 
nm is the interparticle distance. With a particle diame¬ 
ter of 3.7 nm the mass becomes m=5.1.10 -22 kg, leading 
to K ~ 11 N.m -1 consistently with covalent bonds and 
larger than Van der Waals connection probed in soft as¬ 
semblies of nanoparticles mm- 

Beside the evaluation of collective sound velocity, the 
analysis of the coherent acoustic phonons spectrum pro¬ 
vides insights on the ultrafast generation processes. The 
transient optical signal reported for the thicker super¬ 
lattice is clearly different from those obtained earlier for 
films made of noble metals like gold, silver m or copper 
m • The physical characteristics of the photogeneration 
and photodetection processes can be obtained thanks to 
a proper simulation of the transient reflectivity signal 
(AR/R). We performed this simulation for the thicker 
superlattice. We employed the standard transient reflec¬ 
tivity model where in our case AR/2R=Real(£r/r) with 

Q3BG2]: 

Sr/r ~ p + iScf) ~ — 2iko5z (2) 

+ /bfc°nsL dn_sL_ f v ( z ,t)exp(2ik 0 h S Lz)dz. 
\ l ~ n SL) dr l Jo 

where ko = 2tt/X is the probe wave vector in air, 
77 ( 2 , t) is the coherent acoustic phonons strain field prop¬ 
agating perpendicularly to the surface of the sample and 
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Figure 4: (color online) (a) The full (dotted) line curve rep¬ 
resents the photoinduced coherent acoustic phonons spectrum 
with (without) phonon attenuation used to adjust the tran¬ 
sient photoacoustic response shown in Fig. §b). The inset 
shows the emitted coherent acoustic phonon strain (ultrafast 
heating of nanoparticles) with a spatial extension of 42 nm 
used for the simulation. This spatial extension is sketched by 
the red curve in (b) where it is compared to the penetration 
depth of light into bulk gold. 


dhsL/dr) is the photoelastic coefficient. z=0 corresponds 
to the air/superlattice interface. In Eq. [2j the first 
term corresponds to the contribution of the surface dis¬ 
placement (Sz) of the superlattice and the second one 
to the photoelastic contribution (modification of the re¬ 
fractive index induced by the strain field of the coher¬ 
ent acoustic phonons). The strain field has been mod¬ 
eled by applying the standard bipolar shape m with 
77 ( 2 :, t) = —sgn(z — Vst)exp(— \ z — Vst | /£). The 
simulation of the photoinduced strain 77 (z,t) is shown 
as a dotted curve in inset of Fig. ET) where £ is the 
characteristic depth over which the photo-induced strain 
occurs. £ is related to the apparent imaginary part of 
the refractive index k V g™ v at the pump wavelength with 
£=A/ 47 r/^ mp . In order to take into account phonon at¬ 
tenuation (anharmonic damping, defect scattering), we 
have used a cut-off frequency at around 30 GHz (full line 
curves in Fig. |4ja) and inset of Fig. |4ja)) in accordance 
with the experimental FFT (Fig. [2jb)) . The calculation 
takes into account the reflection of the acoustic phonon 
wavepacket at the silicon interface to reproduce the three 
visible acoustic echoes. The best adjusted parameters we 
have obtained concerning the probe optical properties (@ 
2.19 eV) are h p s ™ be = 2.2+0.9i and for the pump (@ 1.59 
eV) /^ mp =1.5. These optical coefficients are plotted 
in Fig. [ljf) and the fitted curve of the time derivative 
d(AR/R)/dt compared to the experimental data is shown 
in Fig. |2|a). The adjustment of the shape of the pho¬ 
toacoustic signal depends also on the photoelastic coeffi¬ 
cient dhsL/dr] (iQj which are unknown for such artificial 


nanoparticles superlattices. However we have tested dif¬ 
ferent scenarios and our adjustment shows that the real 
part of the photoelastic coefficient is dominant over its 
imaginary part. Considering the well known relationship 
dhsL/dr] oc dhsL/dE [10], E being the probe energy, and 
from the energy dependence of optical refractive of the 
NPS refractive index shown in Fig. [ljf ), it appears that 
the derivative of the real part (usl) (in particular @2.19 
eV) is much larger than that of the imaginary part ( ksL )• 
Consequently and in accordance with our simulations, it 
indicates that the real part of the photoelastic coefficient 
may dominate with dhsL/dr] ~ dnsL/dp. 

The results of our simulations confirm that the optical 
properties of the NPS are deeply different than those of 
bulk gold crystal. As a comparison, we have performed 
as well a calculation where we simulate the transient op¬ 
tical signal considering that the photo-induced strain oc¬ 
curs over a depth of 11 nm beneath the surface, as in 
bare gold metal as sketched in Fig. Sb) ( hot electrons 
diffusion is neglected). The simulation is shown in Fig. 
[ 2 |b) (black curve) obtained with an optical refractive in¬ 
dex for gold (@2.19 eV) tlau = 0-25 + 3 i and (@1.59 eV) 
hAu =0.3 + 6i. The discrepancy between this simula¬ 
tion and our experimental result clearly confirms that 
the pump light can deeply penetrate this subwavelength 
plasmonic nanostructure over a distance of about £ ^ 42 
nm, much larger than the optical penetration depth of 
bare gold. This distance over which the conversion of 
light energy into mechanical energy occurs is not com¬ 
patible with transport of hot electrons due the existence 
of insulating contacts in between the nanoparticles as 
already discussed in [25] . As a consequence, the ultra¬ 
fast generation of coherent acoustic phonons in our NPS 
can only be interpreted from an ultrafast heating of the 
NPS assisted by transport of the electromagnetic energy 
through plasmon-polariton propagation perpendicular to 
the surface. 

As a summary, we have demonstrated the gener¬ 
ation and detection of coherent acoustic phonons in 
gold nanoparticles superlattices. The analysis of the 
GHz sound propagation has permitted to estimate the 
nanocontact stiffness. Our results demonstrate that the 
emission of coherent acoustic phonons occurs in the NPS 
over a distance of about 45 nm. The spreading of the 
conversion of light energy into mechanical energy requires 
transportation of the light energy in the subwavelength 
regime. We attribute this effect to an ultrafast heat¬ 
ing assisted by plasmon-polariton propagation perpen¬ 
dicular to the NPS surface. The model of the genera¬ 
tion and detection of coherent acoustic phonons is con¬ 
sistent with this mechanism and the optical parameters 
extracted from our fitting are in agreement with conven¬ 
tional optical characterizations. These results open new 
perspectives for the manipulation of plasmons by coher¬ 
ent phonons or viceversa as recently investigated both 
theoretically and experimentally mm- 
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